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SUMMARY 

The  filtered  mass  density  function  (FMDF)  of  mixture  fraction  and  other  filtered  variables 
used  in  large-eddy  simulation  of  turbulent  combustion  are  studied  using  measurement  data 
obtained  in  turbulent  partially  premixed  methane/air  (Sandia)  flames.  For  SGS  scalar  vari- 
ance small  compared  to  its  mean,  the  FMDF  is  not  far  from  Gaussian  and  the  SGS  scalar  is 
well  mixed.  For  large  SGS  variance,  however,  the  FMDF  becomes  bimodal  and  the  condi- 
tionally filtered  scalar  dissipation  is  bell-shaped,  indicating  the  existence  of  a diffusion  (dis- 
sipation) layer  structure,  which  is  similar  to  the  mixture  fraction  profile  in  the  counter-flow 
model  for  laminar  flamelets.  The  conditionally  filtered  temperature  near  the  stoichiometric 
mixture  fraction  decreases  progressively  with  increasing  SGS  scalar  variance.  Furthermore, 
local  extinction  events  appear  to  occur  mostly  when  the  SGS  scalar  variance  is  large,  sug- 
gesting the  possibility  of  flamelet  extinction.  The  results  suggest  that  the  mixing  regimes 
and  the  associated  mixture  fraction  structure  could  potentially  have  strong  influences  on  the 
combustion  regime  and  extinction/reignition  in  turbulent  nonpremixed  flames. 

TECHNICAL  DISCUSSION 

LES  has  been  recognized  as  a very  promising  approach  to  modeling  turbulent  combustion. 
In  LES  the  subgrid-scale  (SGS)  scalar  mixing  and  the  resulting  instantaneous  distribution 
of  scalar  values  in  each  grid  volume  (i.e. , the  filtered  density  function)  must  be  faithfully 
represented  in  order  to  accurately  predict  the  chemical  reaction  rate.  The  present  study 
investigates  issues  in  using  this  approach,  including  SGS  mixing  and  turbulence-chemistry 
interaction,  which  play  crucial  roles  in  LES  of  nonpremixed  combustion. 

Our  previous  studies  of  the  SGS  mixing  of  a conserved  scalar  in  nonreacting  jets  (Tong 
2001,  Rajagopalan  and  Tong  (2003),  Wang  and  Tong  (2002,  2004),  Wang  et  al.  2004)  have 
shown  that  the  SGS  scalar  at  a fixed  location  has  qualitatively  different  FDF  shapes  and 
structures  depending  on  the  instantaneous  SGS  scalar  variance.  When  the  SGS  variance 
is  small  compared  to  its  mean  value,  the  SGS  scalar  has  close  to  Gaussian  distributions, 
indicating  well  mixed  SGS  scalar  fields.  The  SGS  turbulence  is  in  quasi-equilibrium,  i.e.,  the 
production  of  the  SGS  variance  is  equal  to  or  smaller  than  the  SGS  dissipation  rate.  When  the 
SGS  variance  is  large  compared  to  its  mean  value,  the  conditional  SGS  structure  on  average 
resembles  that  of  a counter-flow  diffusion  flame,  which  is  a model  for  laminar  flamelets.  The 
SGS  scalar  has  bimodal  distributions,  indicating  highly  nonpremixed  SGS  scalar  fields.  In  a 
nonpremixed  flame  this  would  indicate  that  the  fuel  lean  and  rich  regions  of  the  SGS  fields 
are  highly  segregated.  There  is  a sharp  interface  separating  the  two  regions,  across  which 
there  is  a large  scalar  value  jump  (can  be  as  large  as  the  integral-scale  fluctuations).  In 
addition  the  SGS  scalar  is  in  nonequilibrium. 

The  well-mixed  SGS  mixture  fraction  fields  with  close-to-Gaussian  distributions  and  the 
highly  nonpremixed  SGS  fields  with  bimodal  distributions  can  potentially  have  strong  influ- 
ences on  the  flame  structure.  In  this  study  we  investigate  the  SGS  mixing  of  the  mixture 


157 


fraction  in  turbulent  partially  premixed  flames  and  examine  the  effects  of  the  different  SGS 
mixture  distributions  and  structures  on  the  flame  structure.  We  investigate  the  characteris- 
tics of  the  filtered  mass  density  function  (FMDF)  of  the  mixture  fraction 


^l(C;x,£)  = (p{x,t)8(£  - i;x,t))e  = / p(x',t) 8 (f  - £;x\f)G(x  - x')dx',  (1) 


and  the  mixing  term  in  its  transport  equation,  the  conditionally  filtered  scalar  dissipation 
rate, 


dxj  dxj 


(2) 


where  £,  p,  and  D are  the  mixture  fraction,  the  fluid  density  and  the  diffusivity  for  the 
mixture  fraction,  respectively.  The  subscripts  £ and  L denote  conventional  and  Favre  filtered 
variables  respectively.  The  dependences  of  PiL  and  (x|£)f  on  the  different  mixing  process 
will  be  studied.  We  also  analyze  the  conditionally  filtered  temperature  (Tj£)f  and  examine 
conditional  SGS  temperature  samples  to  study  the  effects  of  SGS  mixing  on  flames  structure. 


Experimental  data  and  processing  procedures 

We  use  experimental  data  obtained  in  piloted  turbulent  partially  premixed  methane  with 
a T.3  ratio  of  CH4  to  air  by  volume  (Sandia  flame  D and  E,  (Karpetis  and  Barlow  (2004, 
2004).  Their  measurements  employed  combined  line-imaging  of  Raman  scattering,  Rayleigh 
scattering,  and  laser-induced  CO  fluorescence.  Simultaneous  measurements  of  major  species 
(C02,  02,  CO,  N2,  CH4,  H20,  and  H2),  mixture  fraction  (obtained  from  all  major  species), 
temperature,  and  the  radial  component  of  scalar  dissipation  rate  were  made.  The  mixture 
fraction  is  calculated  using  a variation  of  Bilger’s  definition,  which  has  been  modified  by 
excluding  the  oxygen  terms.  The  length  of  the  imaging  line  is  6.0  mm  with  a resolution  of 
0.2  mm. 

Measurements  of  the  filtered  density  functions  require  spatial  filtering  of  scalar  fields.  In 
this  research,  both  one-dimensional  filtering  will  be  employed.  The  filter  sizes  A employed 
in  this  work  are  3.0  and  6.0  mm. 

Results 

The  FMDF  and  other  SGS  variables  are  analyzed  using  their  conditional  averages.  We 
use  the  Favre  filtered  mixture  fraction, 

(Ol  = ( P§t/{p)i , (3) 


and  the  Favre  SGS  scalar  variance, 

({"2)l  = - (OlM  = - «>i,  (4) 

as  conditioning  variables.  The  SGS  variance  value  normalized  by  its  ensemble  mean  is 
denoted  by  {C/l- 

The  conditional  mixture  fraction  FMDF  for  flame  D for  several  SGS  variance  values 
is  shown  in  Fig.  1.  The  results  for  flame  E (not  shown)  are  similar.  The  Favre  filtered 
mixture  fraction  is  set  to  the  stoichiometric  mixture  fraction,  £s  = 0.35,  to  maximize  the 
probability  of  the  SGS  field  containing  reaction  zones.  For  (£,/2)2  <<  1 the  conditional 
FMDF,  (F^/J  (£)/,,  (£//2)l)>  is  generally  unimodal  and  probably  not  very  far  from  Gaussian. 
This  suggests  that  the  SGS  scalar  is  likely  to  follow  the  Kolmogorov’s  cascade  picture. 
Therefore,  when  the  filter  scale  decreases  the  SGS  fluctuations  become  smaller,  suggesting 
that  the  burden  on  the  SGS  mixing  models  is  reduced. 
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Figure  1:  The  mean  conditional  scalar  FMDF  for  Sandia  flame  D for  (£)/,  = and  several 
(C'2)l  values  (specified  valuex[0.67,  1.49]).  (a)  x/D=15,  A = 6.0mm;  (b)  x/D=15,  A = 
3.0mm.;  (c)  x/D=30,  A = 6.0mm 


For  (C'2)l  close  to  or  larger  than  one,  the  FMDF  becomes  bimodal  and  the  bimodality 
is  stronger  for  larger  SGS  variance.  This  indicates  that  the  SGS  mixture 
fraction  is  highly  nonpremixed.  Therefore,  the  rich  and  lean  mixtures  in  a grid  cell  are 
essentially  segregated.  The  bimodal  FMDF  is  also  similar  to  the  scalar  FDF  for  large  SGS 
variance  observed  in  nonreacting  flows.  The  SGS  scalar  is  not  well  described  by  Kolmogorov’s 
turbulence  cascade  picture. 

In  contrast  to  the  unimodal  FMDF,  the  dependence  of  the  bimodal  FMDF  on  the  filter 
scale  is  much  weaker.  The  FMDF  at  x/D  = 15  for  a filter  size  of  3.0  mm  shown  in  Fig.  1(b) 
has  the  same  characteristics  as  that  for  a 6.0  mm  filter.  Therefore,  large  scalar  value  jumps 
exist  in  the  SGS  scalar  regardless  of  the  filter  scale  as  long  as  the  filter  scale  is  larger  than 
the  scalar  dissipation  length  scale,  which  is  generally  the  case  in  LES.  This  suggests  that 
when  the  SGS  variance  is  large,  the  burden  on  the  mixing  model  to  correctly  model  the  SGS 
mixing  is  not  reduced  when  the  filter  size  is  decreased. 

The  results  show  that  the  statistical  structure  of  the  SGS  mixture  fraction  is  differ- 
ent for  small  and  large  SGS  variance  values.  The  difference  of  the  £ values  between  the 
two  FMDF  peaks  is  often  larger  than  the  equilibrium  reaction  zone  width  in  the  £ space 
for  the  methane  flames  studied  in  this  work  (A « 0.23  (Frank  et  al.  2002).  Therefore, 
such  a mixture  fraction  structure  is  likely  to  support  laminar  flamelets.  This  is  in  contrast 
with  the  well-mixed  SGS  mixture  fraction  for  small  SGS  variance  where  the  turbulence  cas- 
cade is  likely  to  dominate  and  the  dissipation-scale  scalar  fluctuations  are  predicted  by  the 
Kolmogorov-Obukhov-Corrsin  theory.  Therefore,  such  a SGS  mixture  fraction  is  likely  to 
result  in  distributed  reaction  zones. 

The  conditionally  filtered  scalar  dissipation  (xl£)*  (n°t  shown)  also  has  qualitatively 
different  dependences  on  the  mixture  fraction  for  small  and  large  SGS  variance.  For  small 
(C2)l  ^ depends  relatively  weakly  on  £,  except  near  £ = 0 and  1.  This  is  consistent  with 
the  FMDF  being  unimodal  and  not  far  from  Gaussian.  For  large  (C2)l->  (x|£)*  becomes  bell- 
shaped, with  the  maximum  at  the  £ value  where  the  bimodal  FMDF  has  a dip,  indicating  that 
there  is  an  interface  between  the  highly  segregated  SGS  mixture  fraction  regions  and  that 
the  interface  is  sharp  and  is  essentially  a diffusion  layer.  The  FMDF  and  the  conditionally 
filtered  dissipation  results  suggest  that  the  SGS  mixture  fraction  structure  is  similar  to  that 
in  the  counter-flow  model  for  laminar  flamelets. 

The  conditionally  filtered  temperature  ((!T|£)f| (£){,,  {C2)})  for  the  Sandia  flame  D at 
x/D  — 7.5  and  15  are  shown  in  Fig.  2.  For  small  SGS  variance  the  conditionally  filtered 
temperature  is  close  to  the  equilibrium  values.  As  the  SGS  variance  increases,  the  tempera- 
ture begins  to  decreases,  especially  near  the  stoichiometric  mixture  fraction.  The  magnitude 
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of  the  temperature  reduction  is  considerably  smaller  for  x/D  = 7.5  than  for  x/D  — 15  and 
x/D  — 30  (not  shown).  Because  the  jump  in  mixture  fraction  is  large  for  large  SGS  vari- 
ance, the  decreases  in  temperature  under  such  conditional  are  likely  caused  by  straining  or 
extinction  of  laminar  flamelets. 

It  is  important  to  note  that  although  the  conditionally  filtered  temperature  in  Fig.  3 is 
somewhat  similar  to  the  conditional  mean  temperature  obtained  by  using  \s  as  the  condi- 
tioning variablesf?],  conditioning  on  the  SGS  variance  provides  additional  information  than 
conditioning  on  the  dissipation  rate  alone.  This  is  because  the  SGS  variance  enables  condi- 
tional sampling  of  the  structure  of  the  SGS  mixture  fraction,  thereby  allowing  the  effects  of 
the  mixture  fraction  structure  on  the  flame  structure  to  be  studied. 


(b) 


Figure  2:  The  mean  conditionally  filtered  temperature  for  flame  D for  (£)/,  = £ s and  several 
(C'2)l  values,  (a)  x/D=7.5,  A = 3.0mm;  (b)  x/D=15,  A = 6.0mm. 
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